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The 3-loop Beta Function of QCD 
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Using the background field technique, we calculate the 3-loop beta function of lattice SU{N) gauge theories. 
In the pure gluonic case, we present our results, comparing to those recently obtained by Luescher and Weisz. 
We also provide a progress report in the case of QCD with Wilson fermions. 



1. INTRODUCTION 

The relationship between the bare coupling and 
the cutoff is an essential ingredient in lattice cal- 
culations. Usual Monte Carlo simulations require 
the knowledge of this relationship far from the 
critical point, where corrections to asymptotic 
scaling could become relevant. 

In order to check the relevance of these scaling 
corrections, we present here a calculation of the 
first non-universal coefficient of the lattice /3 func- 
tion in pure Yang-Mills theory with Wilson action 

1 , as well as in the presence of Wilson fermions 

2 . This rcnormalization group function can be 
written as 
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where a is the lattice spacing, the lattice bare 
coupling, gr is the renormalized coupling constant 
and /X the subtraction point. It establishes how 
the bare coupling and the cutoff a must simul- 
taneously vary to keep renormalized quantities 
fixed. The first two coefficients in eq. (1) are 
known as they are scheme- independent. Our pur- 
pose is to compute the third coefficient, &2 ■ 

In a recent paper |^ the authors calculated a 
quantity related to this coefficient (in the pure 
gluon case), using a coordinate space method 
for evaluating lattice integrals. In our calculation 
we have used a different technique where super- 
ficially divergent integrands are Taylor expanded 
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in D dimensions. The idea is a generalization 
to two loops of the procedure introduced in . 
Thus, our results for the pure gluon case provide 
an independent check of the calculation in refer- 
ence [||. 

In what follows we will first define the quanti- 
ties we set out to compute, and give a brief de- 
scription of our calculational techniques. We then 
present our results. We discuss the effect of 62 on 
the scaling behaviour of dimensionful quantities. 
We close with some concluding remarks. 

The involved algebra of lattice perturbation 
theory was carried out by means of a computer 
code. The main features of this code were out- 
lined in Q and used to compute the three loop 
perturbative background of the topological sus- 
ceptibility dl, and the three loop lattice free en- 
ergy 0, in pure Yang-Mills on the lattice. For 
the purpose of the present work, this code was 
extended to include form factors, and to handle 
fermionic vertices. 

Aside from the intrinsic interest of the result, 
this calculation serves also as a prototype for nu- 
merous 2-loop calculations involving fermions and 
external momenta, e.g. multiplicative rcnormal- 
ization of fermionic currents. We expect to return 
to this issue in a future report. 

2. GENERAL CONSIDERATIONS 

We performed our calculation using the back- 
ground field method. To facilitate comparison, 
we have adopted the notation of Ref . P| . 
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We want to calculate ly'^^^p), defined through 
the lattice background field two point function 
via: 

j2 r(2.<'^o)(p, ~p)';^l = ^s^''3p'{i - Kp)]/3o' (2) 

/J 

with iy{p) = J2 9o'' ^^^Hp) ■ Knowledge of v^^\p), 
together with its MS counterpart ^'^^''(p) and 
the one loop gluon self energy leads us, through 
use of standard formulae, to Z{g{),^a) {go = 
Z{go, ^a)g{fj.)). Knowledge of the three loop coef- 
ficient of the renormalized /3-function, 62^'^, along 
with the two loop coefficient Z20 in the expansion: 

i 

Z{go, A^a) 1 + X! X! do'^tj MO (3) 
i=i j=i 

gives us the three loop lattice /3 function, through 
the relation: 

b^^bf^- 261Z10 + 60^10 + 260Z20 (4) 

All one loop quantities appearing above are well 

(2) 

known. For the continuum function (p), in 
the presence of Nf fermionic species, we find (in 
the Feynman gauge): 

(p) (167r2)2 =iv2[8p + 577/18 - 6C(3)] (5) 

+Nf[{-3p - 401/36)iV+ (p + 55/12 - 4C(3))/iV] 

{p^ln{^,yp^)). 

Contributions to iy^^\p) come from 35 dia- 
grams in the pure gluonic case (shown in Ref.[^). 
Fermions bring in 18 more two- loop diagrams (the 
ghost diagrams of with fermions replacing 
ghosts) plus 2 one-loop diagrams and 2 diagrams 
with fermion mass counterterms (we work with 
zero renormalized fermion mass) . We manipulate 
these diagrams using a computer code 1^ which 
has now been extended to include fermionic ver- 
tices and external momenta. The code starts with 
a series of operations which can be performed 
on diagrams with arbitrary number/type of ver- 
tices/loops/external legs: Generation of vertices, 
contraction with appropriate reassignment of in- 
dices, symmctrization, reduction of Dirac and 
color indices, simplification of Lorentz structures. 
It is essential that the corresponding algorithms 



be (quasi-)polynomial: The contrary results in as- 
tronomical requirements of CPU and RAM, even 
in the simplest cases. This last constraint turns 
out to be the most exacting one in the construc- 
tion of our algorithms. 

At this point, we apply trigonometry judi- 
ciously to bring our expressions into a "canonical" 
form, separating them also into (a few) superfi- 
cially divergent terms plus the remainder. This 
remainder can be quite sizeable, especially in the 
fermionic case (often thousands of terms) but it 
can be handled with a straightforward Taylor ex- 
pansion in external momenta or, when subdiver- 
gences are present, with appropriate subtractions; 
e.g., in 

(2^)-« J d^k°q ■°k°q-°p /{fk^^kk+p") (6) 

... — ... 2 — ... 2 

the subtraction: l/q+k = l/f+{l/q+k -l/q^) 
leads to a product of one loop integrals with a sin- 
gle Inp, plus a Taylor expandable part. Different 
types of subtractions may be necessary for the full 
evaluation of a diagram; suffice it to say that the 
most cumbersome diagram (a "diamond" with a 
triangle of fermions) splits into ^100 types of ex- 
pressions. 

For the superficially divergent diagrams, we 
have extended beyond one loop the proceduce of 
Ref. Q (analytic continuation to 4— 2e dimen- 
sions, followed by Taylor expansion). This entails 
a lengthy series of additions/subtractions to parts 
of the integrand and to the integration domain (so 

that, e.g., l/k+q = + ., l/q^ = l/g^ + - ■ ■, 
[— 7r/a, Tr/a]"' ~ — . . etc.). Interchanging the 
and limits is more complicated in this 
case. Also, certain intermediate expressions must 
be calculated to 0{e). Proceeding in this way 
we have produced [Q a table of values for all such 
two- loop diagrams, complementing those given in 

i- . . . 

Our numerical integrations were done on finite 
lattices, with subsequent extrapolation to infinity, 
using the form X^m n -^)/-^" lattice 
size). This leads typically to four significant dig- 
its. The coordinate space method of Ref. Q gives 
a much better precision, but cannot be applied as 
is to fermionic diagrams. 
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3. RESULTS 

The contribution Viip) of the ith diagram to 
1^(2) (p) is: 

3 4 
tj.=0 P 

+ C3.1npV(4^)' + C4,^}+0(/) (7) 

where: Cn,i = cf\/N^ + c':^\ + c':^\n'^ for diagrams 

without fermions and CnA — c'n\l^ + ^n\^ ^'^^ 
diagrams with fermions. 

Gauge invariance requires that Co,i = 0; this 
is indeed so, separately for diagrams with and 
without fermions, as we have checked both by nu- 
merical evaluation and algebraically (expressing 
all vertices involved in terms of derivatives of the 
gluon and fermion propagators, and performing 
integration by parts where necessary). Lorentz 
invariance requires: ci,i — 0. Such contribu- 
tions appear in non-tadpole diagrams containing 
one-loop renormalized gluon propagators; again, 
we verified that their sum vanishes, both numer- 
ically and algebraically. The coefficients C2,i of 
the logarithms squared must coincide with those 
of the corresponding continuum diagrams and in- 
deed they do. 

In the pure gluonic case, we compared our re- 
sults to those of Liischer and Weisz. All coeffi- 
cients co,i, ci.i, C2,i, C3,i coincide. As for C4,i, the 

parts cf\ and c'l\ agree for each diagram. The 

coefficients c\ ■ agree for all but 1 diagram (# 5 
of ^), leading to a slightly smaller (<3%) total 

(2) 

for c\'. To this writing, we have not been able 
to find the cause of this discrepancy. 

Thus, comparing with the expression of |^ for 
v^'^\p) without fermions: 

= -iVV(327r4)lnp2 + i/Ar2(3/i28) 
-0.01654462 -fiV20.0074438722 (8) 

we find for the last number in the above: 
0.007230. 

Our results for i/^^^ {p) with Wilson fermions 
are in the final stage of extrapolation/precision 
checks 1^. 



4. CONCLUSIONS 

Eq.(|) for 1/(2) (p) can be used directly to eval- 
uate 6j, as explained. In the absence of fermions, 
we find: 

[IQir'^fh^ = -2143/Ar + 1433. BA^ - 366.2A^3 (9) 

Substituting now in the asymptotic scaling for- 
mula: 

akL = cxp(-l/26o5o)(Mo)"''^/''' • (10) 
{l + gl{hl-b^h^)/2bl + ...) 

the last factor becomes, for 5[/(3): 

(1 + 0.189652) (11) 

Thus, in the region ^ 1 this factor brings about 
a rather sizeable correction to asymptotic scaling. 

Some further computations which can be car- 
ried out without changes in our computer code 
are two loop renormalization of various compos- 
ite operators (fermionic currents, etc.), use of im- 
proved actions and (with some modifications) use 
of staggered fermions. 
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